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Abstract 

Homeostatic and adaptive control mechanisms are essential for keeping organisms structurally and functionally stable. 
Integral feedback is a control theoretic concept which has long been known to keep a controlled variable A robustly (i.e. 
perturbation-independent) at a given set-point A set by feeding the integrated error back into the process that generates A. 
The classical concept of homeostasis as robust regulation within narrow limits is often considered as unsatisfactory and 
even incompatible with many biological systems which show sustained oscillations, such as circadian rhythms and 
oscillatory calcium signaling. Nevertheless, there are many similarities between the biological processes which participate in 
oscillatory mechanisms and classical homeostatic (non-oscillatory) mechanisms. We have investigated whether biological 
oscillators can show robust homeostatic and adaptive behaviors, and this paper is an attempt to extend the homeostatic 
concept to include oscillatory conditions. Based on our previously published kinetic conditions on how to generate 
biochemical models with robust homeostasis we found two properties, which appear to be of general interest concerning 
oscillatory and homeostatic controlled biological systems. The first one is the ability of these oscillators ("oscillatory 
homeostats") to keep the average level of a controlled variable at a defined set-point by involving compensatory changes in 
frequency and/or amplitude. The second property is the ability to keep the period/frequency of the oscillator tuned within a 
certain well-defined range. In this paper we highlight mechanisms that lead to these two properties. The biological 
applications of these findings are discussed using three examples, the homeostatic aspects during oscillatory calcium and 
p53 signaling, and the involvement of circadian rhythms in homeostatic regulation. 
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Introduction 

The biological motivation of this work can be summarized as 
follows: How can homeostatic mechanisms possibly work when 
many or even most of the regulatory processes within a cell are 
based on oscillations? Versions of this question and how oscillatory 
processes participate in homeostatic and adaptive mechanisms 
have been repeatedly asked and discussed [1-5]. Our aim is to 
identify and build homeostatic/ adaptive motifs on a rational basis 
with possible applications within physiology and synthetic biology. 
In this paper we apply control-engineering and kinetic methods 
and show how the classical concept of homeostasis [6,7] is linked 
to oscillatory behavior. We demonstrate how biological oscillators 
can have robust (perturbation-independent) homeostatic/adaptive 
behaviors both with respect to average concentration of a 
regulated variable and with respect to a robust control of the 
oscillator's frequency. By taking three examples, we argue that 
such properties appear closely linked to the controlled period 



lengths of the p53-Mdm2 oscillatory system and circadian rhythms 
[1,8] or to the homeostatic regulation of cytosolic calcium during 
signaling [9]. 

Organisms have developed defending homeostatic mechanisms 
in order to survive changing or stressful conditions by maintaining 
their internal physiologies at an approximately constant level 
[7,10,11]. In this respect, many compounds are tightly regulated 
within certain concentration ranges, because they are essential for 
cellular function, but may lead to dysfunction and diseases when 
their concentrations are outside of their regulated regimes. The 
term "homeostasis" was introduced by Cannon [6,7] to indicate 
that the internal milieu of an organism is regulated within narrow 
limits. The examples Cannon addresses in 1929 [6] are still actual 
research topics, such as the regulations of body temperature, blood 
sugar, blood calcium and blood pH levels [12-15]. Today many 
more homeostatic controlled compounds have been identified, 
including hormones [16], transcription factors and transcription 
factor related compounds [17], cellular ions such as plant nitrate 
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levels [18,19], iron [20], and calcium [21]. The Supplementary 
Material of Ref. [22] contains further examples. 

Because many biochemical processes are oscillatory [1,8,23- 
27], Cannon's definition of homeostasis has been perceived as 
unsatisfactory and various alternative homeostasis concepts have 
been suggested. The term predictive homeostasis [1] has been 
introduced in order to stress the anticipatory homeostatic behavior 
of circadian regulation. Other concepts include allostasis [2,5] to 
focus on the concerted and interwoven nature of the defending 
mechanisms, rheostasis [3] to put emphasis on set-point changes, 
and homeodynamics [4] to stress the nonlinear kinetic behaviors of 
the defending mechanisms as part of an open system. 

The appearance of cybernetics together with system theory [28- 
31] caused an interest to understand homeostasis and biological 
control from the angle of system analysis and control theory [32- 
39] by introducing control-engineering concepts such as integral 
control [22,40-45]. Integral control allows to keep a controlled 
variable (say A) precisely and robustly at a given set-point A set by 
feeding the integrated error back into the process by which A is 
generated [46]. To gain insights how integral control and 
homeostasis may appear in biochemical and physiological 
processes, we started [43] to study two-component negative 
feedback controllers, where one component is the (homeostatic) 
controlled variable A, while the other is the manipulated or 
controller variable £. Each controller consists of the two species A 
and £ and three fluxes, the inflow and outflow to and from E and 
an E -controlled compensatory flux (either inflow or outflow) oiA, 
denoted j CO mp- The compensatory flux compensates for distur- 
bances in the level of A caused by perturbations in other 
uncontrolled inflows /outflows of A. By considering activating or 
inhibitory signaling events from A to E and vice versa, eight basic 
negative feedback configurations (controller motifs, Fig. la) can be 
created [22,47]. Two kinetic requirements leading to integral 
control have so far been identified, one based on a zero-order 
kinetic removal of the manipulated variable E [22,43,48], the 
other on an autocatalytic formation of £ in association with a first- 
order degradation [45]. Fig. lb gives a brief summary of these two 
kinetic approaches by using motif 5 as an example. For details, the 
reader is referred to [22,45]. We feel that this approach provides a 
rational basis to build networks which allow to view the behaviors 
of the individual controllers and to understand emergent 
properties of the overall network. By combining individual 
controller motifs with integral control we previously showed that 
an integrative and dynamic approach to cellular homeostasis is 
possible, which includes storage, excretion and remobilization of 
the controlled variables [19,22,49]. 

In the present study we extend the concept of homeostasis to 
include sustained oscillatory or pulsatile conditions. We show that 
oscillatory homeostats based on the controller motifs in Fig. 1 a can 
maintain robust homeostasis in A. For controllers where E is 
inhibiting the compensatory flux (motifs 2, 4, 6, and 8, Fig. la), the 
frequency can be shown to depend on the level of £ and therefore 
on the applied perturbation strength. In this class of controllers the 
frequency generally increases upon increased perturbation 
strengths; here we use motif 2 as a representative example. For 
the remaining controller motifs the frequency has been found to be 
less dependent upon perturbations. As a representative example 
for this behavior we use motif 5. We further show that robust 
frequency control can be achieved by either using additional 
controllers, which keep the average levels of A and E homeostatic 
regulated, or by using the intrinsic harmonic/quasi-harmonic 
properties of motifs 1 or 5. The biological significance of these 
findings is discussed with respect to the oscillatory signaling of 
cytosolic calcium and p53, as well as the regulating properties of 



circadian rhythms with respect to homeostasis and temperature 
compensation. 

Results 

Kinetic Approach to Implement Integral Control 

We consider the negative feedback motifs in Fig. 1. A general 
condition for integral control can be formulated if the rate 
equation of the manipulated variable E allows for a rearrangement 
in form of two functions g(E) and h(A), and where the integral of 
\/g{E) with respect to E exists and can be written as G{E). Then, 
the set-point in A is determined by the solution of h(A) = 0, i.e. 

E=h(A)g(E) (1) 
Rearranging Eq. 1 and requiring steady state conditions gives: 



g(E) 



= G(E) = h(A) = 0 



(2) 



Eq. 2 has been applied for nonoscillatory steady states with 
g(E) = l by using zero-order kinetic degradation/inhibition of £ 
[22,43] or with g(E) = E by using first-order autocatalytic 
formation and degradation in E [45]. Other functions of g(E) 
may be possible but plausible reaction kinetic mechanisms need to 
be identified. For the sake of simplicity, we consider here that 
integral control is achieved by a zero-order removal of E using 
g(E) = h 

To extend the condition of Eq. 2 to sustained stable and 
marginally stable oscillations, we observe that the integral of the 
periodic reaction rates A and E along a closed orbit c in the 
system's phase space is zero. For E this can be written as: 



E>, = i> Edt= <j) h(A)dt = 0 



(3) 



Dependent on whether A is activating or inhibiting the 
production or removal of E, two expressions for the set-point of 
the oscillatory controller can be derived from Eq. 3. In case A is 
activating (motifs 1, 2, 5, 6) and by assuming first-order kinetics 
with respect to A in the rate equation for E, the set-point of A is 
given by (see Eq. SI in (File SI)) 



<A> C = o A{i)dt= <A> S 



(4) 



where the integral is taken along one (or multiple) closed and 
stable orbit(s) in the system's phase space. With increasing time t, 
the average concentration of A, <A> t , will approach its set-point 
< A >„,, i.e., 



<A>,= 



A(x)d%-* <A> set when /->co (5) 



When A is inhibiting the production or removal of £ (motifs 3, 
4, 7, 8) and assuming (for the sake of simplicity) that the inhibiting 
term has a first-order cooperativity with respect to A with an 
inhibition constant Kf , the following expression is conserved and 
perturbation-independent (see derivation in File SI, Eq. S8): 
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Figure 1 . A basic set of two-component homeostatic controller motifs with two implementations of integral control, (a) Compound A 
is the homeostatic controlled variable and E is the controller or manipulated variable [22]. The motifs fall into two classes termed as inflow and 
outflow controllers, dependent whether their compensatory fluxes j amp add or remove A from the system. In motifs outlined in gray the controller 
compound E inhibits the compensatory flux, while in the other motifs E activates the compensatory flux, (b) middle figure shows a standard control 
engineering flow chart of a negative feedback loop, where the negative feedback results in the subtraction of the concentration of A (blue line) from 
A's set-point (red line) leading to the error (A sa — A). The error feeds into the integral controller (brown box). The controller output (the integrated 
error) is the concentration of E (green line) which regulates the process that creates A. The perturbations which affect the level of A are indicated in 
orange color, (b) left panel shows the structure of negative feedback (outflow) controller 5. The colors correspond to those of the control engineering 
flow chart. For example, the set-point (red) is given by the ratio between removing and synthesis rates of E, while the integral controller (brown) is 
related to the processing kinetics of E, in this case E is removed by zero-order [22,43], (b) right panel shows the same outflow controller (motif 5). 
The only difference is that the integral controller is now represented by a first-order autocatalytic formation (indicated by brown dashed arrow) and a 
first-order removal with respect to E [45]. 
doi:1 0.1 371 /journal.pone.01 07766.g001 



■At) 



= constant 



(6) 



Homeostasis by Oscillatory Controllers 

To illustrate the homeostatic response of the oscillatory 
controllers, we use, as mentioned above, conservative and limit- 
cycle versions of inflow controller motif 2 and outflow controller 
motif 5 as representative examples. These motifs have been 
chosen, because they represent different ways to achieve negative 
feedback and homeostasis of the controlled variable A . In motif 2 
(as in motifs 4, 6, and 8) E inhibits the compensatory flux, while in 
motif 5 (as in motifs 1 , 3, and 7) the compensatory flux is activated 
by E. A limit-cycle version of motif 6 will be used to discuss 
cytosolic Ca 2+ oscillations in terms of a homeostatic mechanism. 



Conservative Oscillatory Controllers. A conservative sys- 
tem is a system for which an energy or Hamiltonian function (H- 
function) can be found and for which the values of H remain 
constant in time. Conservative oscillators show periodic motions 
characterized by that they in phase space do not occur in isolation 
(i.e. they are not limit cycles). For a given //-level h & periodic 
motion (a closed path in phase space) is surrounded by a 
continuum of near-by paths, obtained for neighboring values of 
h [50] . The dynamics of a two-component conservative oscillator 
can be derived from the //-function using the following equations: 



8H 



dH 
M 



(7) 



which are analogous to the Hamilton-Jacobi equations from 
classical mechanics. In general, solutions of these equations are not 
necessarily oscillatory, but here we focus only on the conservative 
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oscillators, which can be derived from the eight controller motifs 
(Fig. la). Dependent on how integral control is implemented, 
some of the conservative oscillators are well-known; they are: the 
harmonic oscillator [51] based on either motifs 1 or 5 (using zero- 
order implementation of integral control; see left panel in Fig. lb), 
the Lotka-Volterra oscillator [45,52,53] also here based on motifs 1 
or 5 (but using the autocatalytic implementation of integral 
control; see right panel in Fig. lb), and Goodwin's oscillator from 
1963 [54] based on motif 2. In the literature the Goodwin 
oscillator comes in two versions, which are both based on motif 2. 
There is a conservative oscillator version from 1963 [54] with two 
components. There is also another version from 1965 with three 
components [55]. The difference between the two versions lies in 
the kinetics of the degradation rates of the oscillators' components. 
In the 1965 three-component version the degradation rates are 
first-order with respect to the degrading species, while in the 
conservative case (1963 version) the degradation rates have zero- 
order kinetics. These kinetic differences change the oscillatory 
behavior of the two systems significantly. To get limit-cycle 
oscillations, it is well-known from the literature [56] that the three- 
dimensional system where the components are degraded by first- 
order kinetics requires a cooperativity of the inhibiting species of 
about 9 or higher. Our results presented here using motif 2 
confirms Goodwin's 1963 results that when components are 
degraded by zero-order kinetics the system can oscillate with a 
cooperativity of 1 with respect to the inhibiting species E. Here we 
also extend Goodwin's results by showing that limit-cycle 
oscillations can be created based on motif 2, but still using a 
cooperativity of 1 with respect to the inhibiting species E (see 
below). 

The following two requirements are needed to get conservative 
oscillations for any motif from Fig. la: (i) integral control has to be 
implemented in the rate equation for E, and (ii) all removal of A 
should either occur by zero-order kinetics with respect to A, or, 
when the removal of A is first (or nth)-order with respect to A, the 
formation of A needs to be a first (or nth)-order autocatalytic 
reaction [45]. When conditions (i) and (ii) are fulfilled, a function 
H(A,E) can be constructed, which describes the dynamics of the 
system analogous to the Hamilton-Jacobi equations from classical 
mechanics, where the form of H depends on the system's kinetics. 
Details on how H is constructed for the various situations is given 
in File SI. 

Fig. 2a shows a reaction kinetic representation of motif 2, which 
is closely related to Goodwin's 1963 oscillator [54]. It was 
Goodwin who first drew attention to the analogy between the 
dynamics of a set of two-component cellular negative feedback 
oscillators and classical mechanics [54]. In this inflow-type of 
controller, increased outflow perturbations (i.e., increased k 2 
values) are compensated by a decreased average amount of £ (i.e., 
< E > , Fig. 2b), thereby neutralizing the increased removal of A 
by use of an increased compensating flux 



jmmp = ki Kf /{Kf + E) 



(8) 



In this way the average level of A, < A > , is kept at its set-point 
Vmax/k* ( see Eq- S5 in the File SI). During the adaptation in 
<A> (when k 2 is changed) the controller's frequency as well as 
the <E> -level are affected. The frequency oj for each of the 
eight conservative oscillators can roughly be estimated by a 
harmonic approximation (see File SI), which in case of motif 2 
(Fig. 2a) is given by (assuming k\ = 0) 



\/kyk 4 -Kf 
Kf + E ss 



(9) 



E S s { = k%Kf jk 2 —Kf) is the steady state of E, which is obtained 
when A = 0 (Fig. 2a). Because the level of < E > is decreasing with 
increasing k 2 values, Eq. 9 indicates, and as shown by the 
computations in Figs. 2b and 2c, that the frequency of the 
oscillator increases with increasing perturbation strengths (k 2 
values) while keeping < A > at its set-point. In fact, the increase in 
frequency upon increased perturbation strengths appears to be a 
general property of oscillatory homeostats, where the manipulated 
variable E inhibits the compensatory flux (for limit-cycle examples, 
see below). 

At high k 2 values, i.e., when the E level becomes lower than 
Kf, the compensatory flux j CO mp approaches its maximum value 
ki. At this stage the homeostatic capacity of the controller is 
reached. Any further increase of k 2 cannot be met by an increased 
compensatory flux and will therefore lead to a breakdown of the 
controller. For discussions about controller breakdowns and 
controller accuracies, see Refs. [22,48]. 

The scheme in Fig. 2d shows outflow controller motif 5, which 
will compensate any inflow perturbations of A (due to changes in 
k\) by increasing the compensatory flux 



Jcomp 



= k v E 



(10) 



When Kfj«A and K^f' «E the oscillator is harmonic and is 
described by a single sine function which oscillates around the set- 
point <A > set = Vf'^/k^ with frequency oj = \fkyk~l and a 
period of 2nj ' \jkyk$. Increased levels in k\ (Figs. 2e and 2f) are 
compensated by increased <E> levels which keep <A> at its 
set-point. Harmonic oscillations can also be obtained for the 
counterpart inflow motif 1 (see Fig. S9 and Eqs. S44-S50 in File 
SI). 

For the harmonic oscillators (motifs 1 or 5) <A> -homeostasis 
is kept by an increase in < E > , which matches precisely the 
increase in the (average) compensatory flux without any need to 
change the frequency. For the other motifs either an increase or a 
decrease in frequency is observed with increasing perturbation 
strengths dependent whether E inhibits or activates the compen- 
satory flux, respectively. 

Limit-Cycle Controllers. The conservative oscillatory con- 
trollers described above can be transformed into limit-cycle 
oscillators by including an additional intermediate, and, as long 
as integral control is present, homeostasis in A is maintained by 
means of Eq. 4 or 6. Fig. 3a gives an example of a limit-cycle 
homeostat using motif 2. Dependent on the rate constants the 
oscillations can show pulsatile /excitable behavior (Fig. 3b). In 
these pulsatile and highly nonlinear oscillations <A> homeosta- 
sis is maintained at the set-point <A > set = V^' x /k 4 , although the 
peak value in A exceeds the set-point by over one order of 
magnitude (Fig. 3b). As already observed for the conservative case, 
an increase in the perturbation strength (i.e., by increasing k 2 ) 
leads to an increase in frequency while homeostasis in <A> is 
preserved (Fig. 3c). 

Similarly, a limit-cycle homeostat of motif 5 can be created 
(Fig. 4a) by including intermediate e and maintaining integral 
control with respect to A. With increasing perturbation strengths 
(k\ values, Fig. 4b), homeostasis in <A> is maintained by 
increasing <E>. Compared to the conservative situation 
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Figure 2. Representation and kinetics of conservative oscillators based on motif 2 and motif 5. (a)-(c) "Goodwin's oscillator" (motif 2). 
Conservative oscillations occur when Kf { «A and E«Kf}"; the latter condition introduces integral feedback and thereby robust homeostasis [22,43]. 
(b) Conservative oscillations in A and E, with A'i =0.0, fc 2 = 1.0, i^ = lxl0- 6 , yt 3 = 6.0, Kf = 0.5, £ 4 = 1.0, K^». = 2.0, tfff = 1 x lO" 6 . Initial 
concentrations: yio = 1.5, Eq = 1.0. At time t = 50.0 kj is changed from 1.0 to 3.0. (c) <A>, <E> , and frequency as a function of the perturbation ki. 
While the frequency increases and <E> decreases with increasing ki, <A> is kept at its set-point V^/k^=2.0. (d)-(f) Harmonic oscillator 
representation of motif 5. Conservative (harmonic) oscillations occur when Kf { «A (or ki =0) and E«Kfj" . (e) Harmonic oscillations in A and E, with 
ki = 1 .0 (the perturbation), k 2 = 0.0, k, = 1 .0, K$ n = 1 x lO" 6 , k 4 = 1 .0, K,f;« = 2.0, and K{Y = 1 x lO" 6 . At time t = 50.0 /c, is changed from 1 .0 to 3.0. 
Initial concentrations: Aq = 1.5, Eq = 1.0. (f) <A>, <E>, and frequency as a function of the perturbation k\. Typical for the harmonic oscillator is the 
constancy of the frequency upon changing k\ values. <E> increases with increasing k\, while <A> is kept at its set-point Pjjj»j/fct=2.0. 
doi:1 0.1 371 /journal.pone.01 07766.g002 
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Figure 3. A limit-cycle model of controller motif 2. (a) Reaction scheme. Rate equations: A=k\ — krA/(K^i+A) + kg-a; 

E=k4-A - V*%,-EI(K%f +£); a=kyKf /{Kf + E) - k 9 a. (b) Homeostatic response of the model for three different perturbations (/c? values). 
For time t between 0 and 50 units, k 2 = 1 .0 x 10 3 , for t between 50 and 1 00 units, k 2 = 2.0 x 10 3 , and for t between 1 00 and 1 50 units, k 2 = 3.0 x 10 3 . In 
the oscillatory case <A> at time / is given as <A>, = (l/t)X J ' a A(t')dt' (ordinate to the right) showing that <A> is under homeostatic control 
despite the fact that A peak values may be over one order of magnitude larger than the set-point, (c) <A>, <E>, and frequency values as a 
function of ki. Simulation time for each data point is 100.0 time units. Note that <A > is kept at <A> set independent of ki. Rate constant values (in 



au): A:, = 1.0, fc 3 = 1.0 x10 s , k 4 



= 1.0, Kf-- 



l.Ox 10" 



K M~- 



1.0, F**=2.0, K% 



l.OxlO" 6 , and fc 9 =2.0. It may further be noted that the 



degradation kinetics with respect to A are no longer zero-order as required in the conservative case (Figs. 2a-c). Initial concentrations in (b): Aq -- 
Eo = 0.3, and ao = 166.17. Initial concentrations in (c) for each data point: Ao = 1.725 x 10~ 6 , Eq = 1.585, and oo = 0.861. 
doi:1 0.1 371 /journal.pone.01 07766.g003 



1.5, 



(Fig. 2f), the frequency now shows both shght decreasing and 
increasing values. However, the overall frequency changes are not 
as large as for motif 2, indicating that similar to the harmonic case, 
the frequency of the motif 5 based oscillator has a certain intrinsic 
frequency compensation on k\ -induced perturbations (Fig. 4c). 

Robust Frequency Control and Quenching of Oscillations 

In this section we present for the first time biochemical models 
that can show robust (perturbation-independent) frequency 
control. There are several biological oscillators where the 
frequency/period is under homeostatic regulation. Probably the 
best known example is the temperature compensation of the 
circadian period, i.e. these rhythms show an approximately 
constant period length of about 24 h at different but constant 
temperatures [57]. Temperature compensation is also observed in 
certain ultradian rhythms [58,59]. Another biological oscillator 
with a fairly constant period is the p53-Mdm2 system [60], where 
the number of oscillations may indicate the strength of the DNA 
damage in the cell [61]. 

We show two ways how robust frequency control can be 
achieved. One is due to the presence of quasi-harmonic kinetics, 
i.e. the system, although still being a limit-cycle oscillator, behaves 
more like a harmonic oscillator. On basis of experimental results, 
we believe that the p53-Mdm2 system falls into this category (see 
discussion below). In the other approach, frequency homeostasis is 



obtained by regulating E itself by additional inflow/outflow 
controllers I\ J2 ■ This approach leads to many possible ways how 
I1J2 can interact with the central negative feedback A-E loop/ 
oscillator and several ways are illustrated using motif 2 and motif 
5. Such an approach may apply to the period homeostasis of 
circadian rhythms (see discussion below). 

Robust Frequency Control by Quasi-Harmonic 
Kinetics. We consider now the case when the intermediate 
that has been implemented to obtain limit-cycle behavior 
(compounds a or e in Figs. 3a or 4a) obeys approximately the 
steady-state assumption, i.e., axO or e«0. We term the 
oscillators' resulting behavior as quasi-conservative, because these 
systems still have a limit-cycle, but behave also as a conservative 
system. An interesting case occurs when the system is quasi- 
harmonic, i.e. when motifs 1 or 5 are used. In this case the limit- 
cycle oscillations and the frequency can approximately be 
described by a harmonic oscillator, i.e., a single sine function. 
This is illustrated in Fig. 5 where an increased ks value is applied 
to the scheme of Fig. 4a (which leads to e«0). Fig. 5a shows the 
oscillations for three different perturbations [k\ values). The 
oscillations in A show a practically perfect overlay with a single 
sine function, outlined in black for k\ = 1 .0. When k\ is increased 
the oscillations (outlined in blue) undergo a phase shift and an 
increase in amplitude, but the frequency stays constant at the value 
of the (quasi) harmonic oscillator. For high k\ values the 
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Figure 4. A limit-cycle model of controller motif 5. (a). Rate equations: A = k\—krE-A/(K^ l +A) — kyA; e = k 4 -A — k^-e; 
E = k 5 e-V^-E/(K^" +E). (b) Homeostatic behavior in <A> illustrated by three different perturbations (k\ values). At time t= 500.0 k\ is 
changed from 4.0 to 10.0, and at t= 1000.0 k\ is changed from 10.0 to 20.0 (indicated by solid arrows). The set-point of <A> is given as 
FfS/fct=2.0. Rate constant values: k\ is variable, /c 2 = 1.0, iS^ = 0.1, /c 3 =0.0, /c 4 = 0.5, /c 5 = 0.2, V%«=1.0, and X§ B = 1.0xlO -6 . Initial 
concentrations: A$ = 1.9964 x 10~ 2 , eg = 8.0983, and Eq = 12.0258. (c) <A>, <E>, and frequency values as a function of k\ showing that <A> is 
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.Eo = 18.8155. Simulation time for each data point is 10000.0 time units. 
doi:1 0.1 371 /journal.pone.01 07766.g004 



^4-amplitude of the oscillator becomes saturated, which is a 
secondary effect of the oscillator's homeostatic property. Due to 
symmetry reasons and because the oscillator is locked on to the 
harmonic frequency, the value of A cannot exceed beyond twice 
the level of its set-point, which in this case has been set to 12.5 
(Figs. 5a and 5b). As in the harmonic case (Fig. 2f), <E> 
increases with increasing k\ (Fig. 5b). Fig. 5c shows the approach 
to the limit-cycle (outlined in black). When kj increases further 
and the steady state approximation for e becomes better and 
better, the limit-cycle disappears and the system becomes purely 
harmonic. 

Quenching of Oscillations in Quasi-Conservative 
Systems. A requirement to obtain conservative oscillations 
and an oscillatory promoting condition for limit cycle oscillations 
is the presence of zero-order degradation in A. Changing the zero- 
order degradation in A may lead to the loss of oscillations. For 
example, in quasi-conservative systems the oscillations can be 
effectively quenched by either adding a first-order removal term 
with respect to A (with rate constant k^, Fig. 4a) or by replacing 
the zero-order kinetics degradation in A (using kj, Kfj) by first- 
order kinetics with respect to A, or by increasing Kj^. Fig. 5d 
illustrates the suppression of the quasi-harmonic oscillations by 
adding a first-order removal with respect to A. In contrast, when 
an oscillatory system does not show quasi-conservative kinetics, 
addition of a first-order removal with respect to A does not 
necessarily abolish the oscillations. A detailed parameter analysis 



showing how the value of ks affects the period of the oscillations 
and how first-order degradation in A affects the size of the 
parameter space in which sustained oscillations are found is given 
in (Figs. S10 and SI 1 in File SI). 

Robust Frequency Homeostasis by Control of 
<E> . When considering the relationship between <E> and 
the frequency, as for example shown in Fig. 3c, we wondered 
whether it would be possible to design an oscillator with a robust 
frequency homeostasis by using an additional control of < E > . 
For this purpose, two extra controllers I\ and I2 with their own set- 
points for <E> are introduced. Note, that the integral control for 
< A > by £ is still operative and has its own defined set-point. In 
the following we show three examples of robust frequency control 
using motifs 2 and 5. Two of the examples illustrate different 
feedback arrangements of I\ and 1% using motif 2. An example 
using still another arrangement using motif 2 is described in File 
SI (Figs. S12-S14). 

In Fig. 6a a set-up for robust frequency homeostasis is shown by 
using a limit-cycle oscillator based on controller motif 5. The set- 
points for <E>, given by the rate equations for I\ and 1%, are 
< E >set=h/V% ax and <E >£,= V^/k-j. Fig. 6b shows the 
results for a set of calculations when k\ varies from 1 to 20 au. In 
these calculations it was assumed that the I\ and h controllers 
have the same set-point of 20.0 au. In the absence of controllers I\ 
and 1%, the frequency varies as indicated in Fig. 4c, which in 
Fig. 6b is shown as gray dots. When I\ and I2 controllers are both 
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Figure 5. Quasi-harmonic behavior of motif 5 oscillator (Fig. 4a). For time t < 300, a perfect overlay between the numerical calculation of A(t) 
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adjusted to give a closely matching overlay. Other rate constant values (numerical calculations): /r7 = 5.0x 10~ 2 , £3=0. 0, AT^ = 1.0x 10 -6 , fci = 0.8, 
/c 5 = 20.0, K^=10.0, and K%* = 1.0x 10~ 6 . Initial concentrations: A 0 = 12.4290, e 0 = 0.4952, and E a = 1.0139 x lO" 4 . At times Z = 300 and ? = 600 
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the set-point, (b) <A> , A imp i, <E>, and frequency as a function of k\. Simulation time for each data point is 1000.0 time units, (c) Demonstration 
of limit-cycle behavior of the quasi-harmonic oscillations. Same initial conditions as in (a) with k\ = 1 .0, and eo = 0.4952. (d) Same system as in (a), but 
at times t = 200 and / = 500 (solid arrows) ki is changed and kept to 0.1 . The oscillations are efficiently quenched, but A remains under homeostatic 
control. 

doi:1 0.1 371 /journal.pone.01 07766.g005 
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active, <E> shows robust homeostasis at 20.0 (Fig. 6b) and the 
frequency is practically constant (black dots). Fig. 6c shows the 
response when controller I\ has been "knocked out". While in this 
case the < A > values are still under homeostatic control, <E> 
approaches its set-point (defined by <E>l 2 el ) only at high k\ 
values, but without a control of the frequency. When controller I2 
is knocked-out (Fig. 6d), control of <E> and frequency 
homeostasis is observed. Due to the absence of controller I2, 
homeostasis in <E> and in the frequency is lost for higher k\ 
values. The role of I2 in this type of regulator is to diminish/ 
suppress the inflow to A by k\ , such that controller I\ can supply 
the necessary amount of A in order to keep <E> and the 
frequency under homeostatic control. This mechanism is illustrat- 
ed in Fig. 6e by a "static" work mode of h, where the 
concentration of I2 is kept constant. In this case the k\ -region of 
frequency homeostasis increases with increasing but constant 
concentrations of I2 (Fig. 6f). 

A corresponding approach to achieve robust frequency homeo- 
stasis by using motif 2 is shown in Fig. 7a. The set-up differs from 
that used for motif 5 (Fig. 6a) by allowing that I\ and I2 act upon a 
and upstreams of A. For the sake of simplicity, both controllers are 
assumed to have set-points at 20.0 au. Note that in this version of 
the motif 2 oscillator, the removal of A is now purely first-order 
with respect to A (using only ki). Because motif 2 has been the 
core for many circadian rhythm models, we will below discuss 
implications of robust frequency control with respect to properties 
of circadian rhythms. In this context we note that the region 
outlined in gray in Fig. 7a shows the part of the oscillator where 
rate constants have no influence on the frequency, i.e. the 
sensitivity coefficients difrequency) / 8ki are zero. 

Fig. 7b shows the homeostatic behavior in frequency (black 
dots) in comparison with the uncontrolled oscillator (gray dots). In 
the controlled case, both <A> and <E> are under homeostatic 
regulation with set-points of 2.0 au and 20.0 au, respectively. To 
elucidate the effect of the added controllers I\ and I2, we removed 
them one by one (knocking them out). In Figs. 7c and 7d 
controllers I\ and I2 have been removed, respectively. When 
outflow controller I\ is not operative, the system is not able to 
remove sufficient a at low ^2 values. In this case <E> levels are 
high and unregulated at low A^'s and showing an increase in 
frequency. Only at sufficiently high ^2 values controller I2 is able 
to compensate for the decreased levels in <E> . The situation is 
reversed in Fig. 7d, when controller I2 is not operative. At low ^2 
values controller I\ can remove excess of E by diminishing the 
level of a and keeping < E > at its set-point. However, the <E> 
regulation breaks down at high values of k2, because no additional 
supply for E via a can now be provided. In this way controllers I\ I 
I2 act as an antagonistic pair of outflow/inflow controllers, 
respectively. Note that the by E controlled level of <A> (with 
set-point of 2.0 au) is kept at its set-point independently whether 
<E> is regulated by I\ II2 or not. Fig. 7e shows the oscillations 
when both I\ and I2 are operative (Fig. 7b, black dots) and ^2 
being changed from 3.0 to 8.0 at t = 250.0 units (indicated by 
arrow). The level of <E> is controlled to its set-point (20.0), 
while the amplitude of A has increased with the increase of &2. For 
each spike (after steady state has been established) the average 
amount of A is the same and independent of the value of fe, 
leading to the same frequency and homeostasis in < A > . 

Oscillator with Two Homeostatic Frequency Domains 

In the I\ and /2-controlled oscillators described above the set- 
point of <E> will determine the frequency. Fig. 8a shows an 
example of a motif-2-based homeostat, where I\ and h feed back 
to A and a, respectively. For an example where I\ and I2 feed back 



to A only, see Fig. S12 in File SI. In the calculations of Fig. 8, 
different set-points for <E> by controllers I\ and I2 have been 
chosen. As a result, dependent whether the perturbation strength 
(value of k2) is high or low, the oscillator shifts between two 
different homeostatic controlled frequency regimes separated by a 
transition zone (Fig. 8b). Fig. 8c shows the oscillations, <A> and 
<E> values and the frequency switch when k2 is changed from 
3.0 to 8.0. 

Discussion 

Classifications of Biochemical Oscillators and Influence of 
Positive Feedback 

There has been several approaches how chemical and 
biochemical oscillators can be understood and classified [62-66]. 
The controller motifs shown in Fig. la can be considered as a 
basic set of negative feedback oscillators. For example, the Lofka- 
Volterra oscillator can be viewed as a negative feedback oscillator 
based on motifs 1 or 5, but where integral control is implemented 
in terms of autocatalysis [45] and where the controlled variable A 
is formed by autocatalysis and degraded by a first-order process 
with respect to A. The same motif can show harmonic oscillations, 
when integral control and removal of the controlled variable is 
incorporated by means of zero-order kinetics. Two additional 
oscillator types based on the same motif can be created by 
implementing mixed autocatalytic/zero-order kinetics for integral 
control and for the generation/degradation of the controlled 
variable ('Text SI'). The other motifs can be extended in a similar 
way, giving rise to 32 basic (mosdy unexplored) oscillator types. 
This type of classification supplements the one given earlier by 
Franck, where the eight negative feedback loops where combined 
with their positive counterparts to create what Franck termed 
antagonistic feedback [63] An often discussed question is the role 
positive feedback, or autocatalysis, may play in biological 
oscillators. Using a Monte-Carlo approach Tsai et al. [26] studied 
the robustness and frequency responses of oscillators with only 
negative feedback loops and oscillators with a combined positive- 
plus-negative feedback design. The authors concluded that the 
combination of a negative and a positive feedback is the best 
option for having robust and tunable oscillations. In particular, the 
positive loop appears necessary to make the oscillator tunable at a 
constant amplitude. We here have shown how homeostasis and 
tunable oscillators may be achieved without any positive feedback 
(but generally associated with a changing amplitude). To put our 
results in relation to those from Tsai et al. [26], we wondered, 
triggered by the comments from a reviewer, how an oscillator with 
an autocatalytic-based integral controller might behave in 
comparison. For this purpose we used controller motif 2 (Fig. 9a), 
analogous to the scheme shown in Fig. 3a. Interestingly, and in 
agreement with the findings by Tsai et al. [26], the autocatalytic 
step resulted now in relaxation-type of oscillations. As expected, 
the frequency of the oscillator increases with increasing perturba- 
tion strengths &2, and <A> is under homeostatic control 
(Fig. 9b). However, as indicated by the results of Tsai et al. the 
oscillator's amplitude has now become independent of ^2! These 
results show that Franck's original concept of antagonistic 
feedback, i.e. combining positive and negative feedback loops in 
various ways [63] appear to be of relevance for many biological 
oscillators [26]. 

Homeostatic Regulation under Oscillatory Conditions 

In his definition of homeostasis Cannon introduced the term 
homeo instead of homo to indicate that certain variations in the 
concentrations of the homeostatic controlled species are still 
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Figure 6. Oscillator based on motif 5 with robust frequency control, (a) Reaction scheme. Rate equations: 

A=ki K?/(Kp +I 2 ) + k g -h ~krE-AI(Ki, + A)-kyA; e = k 4 A-k 5 e; E=kye- V^E/{K^ + E); h =k 6 -E-V% ax -I 1 /(K*+I 1 ); h = k r E- 
^max'k/(J^M + h)- (b) Demonstration of robust frequency control. <A>, <E>, and frequency are shown as functions of k\. Rate constants: 
Jfci =1.0, k 2 =1.0, A";* = 0.1, fc 3 = 0.0, fc, = 0.5, /c 5 =0.2, Pjjjj" = 1.0, K E g- = 1.0 x lO" 6 , A' 6 = 20.0, (^ = 1.0, = 1.0 x lO" 6 , fc 7 = 1.0, ^=20.0, and 
X^ = 1.0x 10~ 6 . Set-points for E by controllers /, and h are given as <£>{;,, = A 6 /f / ,f; av = 20.0 and <£>^,= Vj* a /ky=20.0, respectively. Initial 
concentrations for each data point (black dots): A 0 = 0.7638, £ 0 = 18.8155, e 0 = 1.6887, 7i, 0 = 1.6695 x 10 3 , and 7 2 ,o = 2.7657 x 10 2 . Gray dots show 
the frequency as a function of k\ without control by l\ and li. (c) System as in (b), but controller l\ not present, (d) System as in (b), but controller h 
not present, (e) Reaction scheme of oscillator, but with a constant h concentration. Rate constants otherwise as in (b). (f) Frequency as a function of 
k\ for the system described in (e) using different constant I 2 concentrations (indicated within the graph). The homeostatic region of the frequency 
increases with increasing I 2 concentrations. 
doi:1 0.1 371 /journal.pone.01 07766.g006 



allowed, but within certain limits [6] . As typical examples, Cannon 
mentions the variations of body temperature, variations in blood 
sugar, blood calcium, and blood pH levels [6]. We have shown 
that the concept of homeostasis can be extended to oscillatory 
conditions and that the term set-point still can be given a precise 
meaning, even when peak values of the controlled variable may 
exceed the set-point by over one order of magnitude (Figs. 3 and 
7). In these cases the set-point relates to the mean value of the 
oscillatory species, < A > . Many compounds are known to be 
under a tight homeostatic regulation to avoid cellular dysfunction, 
such as is the case for cytosolic calcium. There is no particular 



reason to assume that protective homeostatic mechanisms should 
cease to exist once a compound becomes oscillatory and functions, 
as in case of calcium, as a signaling device. Allowing a species (such 
as cytosolic calcium) to oscillate while defending the mean value of 
these oscillations makes it possible to relay signaling without 
exposing the cell to long term overload. In the following we discuss 
three examples where oscillatory homeostats appear to be 
involved: in the homeostatic regulation of calcium and p53 during 
oscillations/ signaling, and in the homeostatic function and period 
regulation of circadian rhythms. 
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Set-points for E by controllers 1\ and I 2 are given as <£>'),,= V^ un jk\\ =20.0 and <£>';,, =/ci 4 / K'j av = 19.998, respectively. The set-point 
<-£■>'),, of the outflow controller Ii has been set slightly higher than <£>';,, for the inflow controller I 2 to avoid integral windup and that the 
controllers work "against" each other [22]. Initial concentrations for each data point (black dots): Ao = 50.4903, Eq =23.9425, ao = 3.2629, 
/i,o = 8.29 55 x 10 3 , and /?,o = 57.8533. Gray dots show the frequency as a function of k 2 for the uncontrolled case, i.e., in the absence of controllers I\ 
and I 2 . (c) System as in (b), but controller I\ is "knocked out" by setting k\\ and I\$ to zero. Homeostasis occurs only at high k 2 values when 
controller I 2 is active, (d) System as in (b), but inflow controller I 2 is inactivated by setting ku and /2.0 to zero. Frequency homeostasis is observed for 
low k 2 when controller I\ is active. At high k 2 values the frequency homeostasis breaks down, because controller I 2 is not present to compensate the 
increased outflow of A, which leads to low <E> values, (e) Oscillations of system in (b) illustrating frequency homeostasis. At time / = 250 (solid 
arrow) k 2 is changed from 3.0 to 8.0. Initial concentrations: A 0 = 21.3167, £0 = 31.7283, a 0 = 0.1237, / L o = 1.3473x 10 4 , and / 2 . 0 = 5.0919 x 10 2 . 
doi:1 0.1 371 /journal.pone.01 07766.g007 
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Figure 8. Oscillator based on motif 2 with robust frequency control but alternative feedback regulation by /] and /> (a) Reaction 

scheme. Rate equations: A = kra-krA/yC& l +A)-k°-A-Ii/(K& 2 +A), E = k 4 - A - V%£.- E / (K^' + E); a = kyKf /(Kf + E)-k 9 -a + k' g I 2 ; h = 

k n -E- V^-hKK'h+hy, t 2 =k iA -E- V^-h/(K^ + / 2 ). (b) Using different set-points = V^/ku =10.0 and <E> , £ l = k u /V% llx = 5.0, 

the frequency (solid dots) can switch between two homeostatic frequency regimes, dependent whether k 2 is low or high. The two regimes are 
separated by a transition zone. Rate constants: /fc 2 = 1.0, fc 3 =20.0, /c 4 = 0.1, V^. = \.5, Kjf = 1.0xl0- 6 , fcn = 1.0, V% ax = 10.Q, x£ = 1.0x 10 -6 , 
/c 14 = 5.0, V% ax =l.O, Z^ = 1.0xl0- 6 , Kf = 1.0, /^=^ = 1.0xl0- 2 , = K^ 2 = 1.0 x IQ- 6 . Initial concentrations: yl 0 = 0.6677, £ 0 = 1.0536, 
a 0 = 2.5828 x 10~ 2 , / L0 = 1.1614 x 10 3 , and / 2 ,o = 7.5008 x 10 2 . (c) Oscillations of system in (b) illustrating frequency switch. At time ( = 500 (solid 
arrow) k 2 is changed from 2.0 to 8.0. Initial concentrations: = 19.7178, E 0 = 0.6272, a 0 = 0.4 1 78, /,, 0 = 1.3 6 96 x 10 2 , and 7 2 ,o = 12.9828. 
doi:1 0.1 371 /journal.pone.01 07766.g008 



Calcium Signaling. Cytosolic calcium (Ca~ + ) levels are 
under homeostatic control to concentrations at about 100 nM 
while extracellular levels are in the order of 1 mM. High Ca 2+ 
concentrations are also found in the endoplasmatic reticulum (ER) 
and in mitochondria (between 0.1-10 mM), which act as calcium 
stores. To keep cytosolic Ca 2+ concentrations at such a low level 
Ca 2+ is actively pumped out from the cytoplasm into the 
extracellular space and into organelles by means of various Ca 2+ 
ATPases located in the plasma membrane (PMCA pumps) and in 
organelle membranes [21,67]. Dysfunction of these pumps leads to 
a variety of diseases including cancer, hypertension, cardiac 
problems, and neurodegeneration [68-70] . During Ca 2+ signaling 
[71,72] cytosolic Ca 2+ levels show oscillations [73-75] but 
signaling can also occur as individual sparks or spikes [76]. Ca 2+ 
oscillations have been found to occur in many cell types and differ 
considerably in their shapes and time scales with peak levels up to 
one order of magnitude higher than resting levels. Similar to the 
behavior of stimulated (perturbed) oscillatory homeostats as for 
example shown in Fig. 3b, Ca 2+ oscillations have been found to 
increase their frequency upon increased stimulation of cells [73— 



75]. The frequency modulation of Ca oscillations [77] is 
considered to be an important property for controlling biological 
processes [75]. The tight homeostatic regulation of cytosolic 
calcium combined with its oscillatory signaling suggests that 
oscillatory homeostats appear to be operative also under signaling 
conditions. 

Although a variety of mathematical models have been suggested 
to describe Ca oscillations [78-84], none of them have so far 
included an explicit homeostatic regulation of cytosolic Ca 2+ . 
Fig. 10a shows how Ca + oscillations can be obtained based on an 
outflow homeostatic controller, which removes excess and toxic 
amounts of cytosolic Ca + . The model considers a stationary 
situation of an activated cell, where a Ca + channel is activated by 
an external signal leading to the inflow of Ca + into the cytosol. 
The increased Ca 2+ levels in the cytosol induce an additional 
inflow of Ca + from the internal Ca + store, a mechanism termed 
"Calcium-Induced Calcium Release" (CICR) [85]. Both inflows 
are lumped together and described by rate constant k\ . The CICR 
flux is maintained by pumping cytosolic Ca 2+ into the ER and 
keeping the Ca 2+ load in the ER high. It should be mentioned that 
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time (au) 

Figure 9. A limit-cycle model of controller motif 2 using autocatalysis as an integral controller, (a) Reaction scheme. Rate equations: 
A=k l -k 2 -A/(K^+A) + k 9 -a; E = k A A-E-k b E; a = kyKf /(Kf + E)-k 9 -a. (b) Homeostatic response of the model for three different 
perturbations (kz values). For time / between 0 and 250 units, /c2 = 0.5, for / between 250 and 500 units, /c2 = 2.0, and for t between 500 and 750 
units, = 5.0. <A > at time t is defined as in Fig. 3. (c) <A>, <E>, and frequency values as a function of ki. Simulation time for each data point is 
2000.0 time units. Note that <A> is kept at A sa =\S.Q (solid black line) independent of k 2 . Rate constant values (in au): k\ =0.0, kj =20.0, k 4 = 0.1, 
Kf = 1.0, Kf { = 1.0 xl0- } ,k 6 = l.5, and /c 9 = 30.0. Initial concentrations in (b): A 0 = 22.09, E 0 = 1.71 x 10 10 , and a Q = 4.0 x 10-". Initial concentrations 
in (c) for each data point are the same as in (b). 
doi:1 0.1 371 /journal.pone.01 07766.g009 



the cause of the Ca + entry across the plasma membrane into the 
cytosol is not fully understood and different views have been 
expressed how this can occur [86,87]. 

For the sake of simplicity, the Ca 2+ concentration in the ER is 
considered to be constant and only the pumping of Ca 2+ from the 
cytosol into the extracellular space is taken into account without an 
increased cooperativity (Hill-function) with respect to the Ca + 
concentration. Fig. 10b shows the oscillations of cytosolic Ca 2+ 
and the homeostat's performance at different inflow rates k\-Ca 2 e ^ t 
into the cytosol, which can reflect different external Ca 2+ 
concentrations and/or different activation levels of the cell. As 
observed experimentally [74] the period of the oscillations 
decreases with increased external Ca 2+ concentration or with an 
increased stimulation of the cell. As shown by < Cfl^+ > t in 
Fig. 10b and by total <Ca 2 c ^ t > in Fig. 10c, on average, robust 
Ca + homeostasis is preserved at varying Ca + inflow rates. In the 
absence of oscillations the Ca + concentration is still kept at its 
homeostatic set-point (Fig. 10b). 

Why Ca 2+ oscillations? A non-oscillatory signaling mechanism 
by cytosolic Ca + would clearly be limited, because a homeostatic 
regulation of cytosolic Ca 2+ would not allow varying Ca 2+ levels as 
a function of external stimulation strengths. On the other hand, a 
frequency-based signaling due to an oscillatory Ca + -homeostat 
would overcome these limitations, because homeostasis is still 
maintained. This has been a brief outiine on how Ca + oscillations 
may be understood on basis of oscillatory homeostasis. More 



detailed studies will be needed, for example by including the 
homeostatic aspect in existing models in order to investigate in 
more detail the implications oscillatory homeostats have on the 
regulatory role of Ca 2+ . 

p53 Signaling. p53 is a transcription factor with tumor 
suppressor properties. In more than half of all human tumors p53 
is mutated and in almost all tumors p53 regulation is not 
functional [88]. In the presence of DNA damage and other 
abnormalities p53 initiates the removal of damaged cells by 
apoptosis. A central negative feedback component in p53 
regulation is Mdm2, an ubiquitin E3 ligase, which leads to the 
proteasomal degradation of p53 and other tumor suppressors [89]. 
In the presence of DNA damage, p53 is upregulated by several 
mechanisms [90-92], and both p53 and Mdm2 have been found 
to oscillate [60] . An interesting feature of these oscillations is that 
their amplitude is highly variable, while their frequency is fairly 
constant [60] . The mean height of the oscillations was found to be 
constant [61]. It was also found that with an increased strength of 
DNA damaging radiation the number of cells with increased p53 
cycles increased statistically [61]. Jolma et al. [51] used the basic 
negative feedback motif 5 (where A is p53 and E is Mdm2) and 
found that the influence of noise on the harmonic properties of the 
oscillations was able to describe the variable amplitudes and the 
approximately constancy of the period. Fourier analysis of the 
experimental data indeed showed that the p53-Mdm2 oscillations 
have a major harmonic component [93] supporting a quasi- 
harmonic character of the p53-Mdm2 oscillations. For such 
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Figure 10. A homeostatic model of cytosolic Ca 2+ oscillations. The model considers a stimulated non-excitable cell under stationary 
conditions using an extended version of outflow controller motif 6, where E is the controller molecule. Intermediate e has been included to get limit- 
cycle oscillations. Rate constant k\ describes the total inflow of Ca 2+ from the ER and from the extracellular space into the cytosol and reflects the 
strength of the stimulation. For the sake of simplicity the external Ca 2+ concentration (Ca 2 J) is considered to be constant (Ca 2 + =1.0). Ca 2 ^ 
denotes cytosolic Ca 2+ and its concentration, (a) Reaction scheme. Rate equations: Ca^ t =k\—k2-KvCc^* t /((K^+Ca^yfKi+E)); 
e = ki-k 4 -e-Ca 2 c +/(K$ + e); E=k 5 e-k f> E. Rate constants: k\, variable; /t 2 = 500; fc 3 =2.0; £4 = 1.0; K$ = 1.0 x 10- 6 ; fc 5 =fc 5 = 1.0; K,=0.l. The 
homeostat's set-point for Ca 2 + is given by £3/^4 = 2.0. (b) Ca 2 + oscillations and average cytosolic Ca 2+ concentration, <Ca 2 + >,, at different 
stimulations and as a function of time /. Initial concentrations: Ca 2 + 0 = 1.772, eg = 2.908 x 10~ 3 ; Eg = 1.643. The quenching of oscillations at low /q is 
due to an increased Km value, (c) Period length and average cytosolic Ca 2+ concentration (<Ca 2 + >) calculated after 2000 time units for different 
stimulation strengths (k\ values). Same rate constants as in (b) with Km = Q.0\. Initial concentrations for each calculated data point: 
Ca 2 + 0 = 6.126x 10- 2 , e 0 = 30.693; £ 0 = 28.806. 
doi:10.1371/journal.pone.0107766.g010 



harmonic or quasi-harmonic oscillations our results (Figs. 2f and 
5b) indicate that p53 is homeostatic regulated both in average 
concentration and in period length to allow to expose the system 
probably to an optimum amount of p53 during each cycle. 
Because the number of p53 cycles appear positively correlated 
with an increased exposure of damaging radiation, the total 
amount of released p53 may be related to a repair mechanism. A 
support along these lines comes from a recent study, which 
indicates that p53 oscillations lead to the recovery of DNA- 
damaged cells, while p53 levels kept at their peak value lead to 
senescence and to a permanent cell cycle arrest [94] . Thus, like for 
cytosolic Ca" + , elevated and oscillatory p53 levels seem to remain 
under homeostatic control in order to mediate signaling events 
and information which appear to be encoded in the oscillations. 

Homeostasis of the Circadian Period. Circadian rhythms 
play an important role in the daily and seasonal adaptation of 



organisms to their environment and act as physiological clocks 
[8,95,96]. Functioning as clocks, their period is under homeostatic 
regulation towards a variety of environmental influences, such as 
changing temperature ("temperature compensation") or food 
supply ("nutritional compensation"). Circadian rhythms partici- 
pate in the homeostatic control of a variety of physiological 
variables, such as body temperature, potassium content, hormone 
levels, as well as sleep [1,8,95,96]. As an example, potassium 
homeostasis in our bodies is under a circadian control, where 
potassium ion is daily excreted with peak values at the middle of 
the day [1]. 

One of the questions still under discussion is how the circadian 
period P is kept under homeostatic control as for example seen in 
temperature compensation. In the antagonistic balance approach 
[97] the variation of the period P with respect to temperature T, 
expressed as d In P/d In T, is given as the sum of the control 
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coefficients [36] C£ = 3 InP/d ln/tj multiplied with the RT- 
scaled activation energies E, (R is the gas constant): 

P _dlnP_^/ain/>\ fd\nk i \_ sr P ( E t \ 
C T~ dlnT ~ ^{sinkJXdlnT) ~ 2f Ck i\RTj 

The sum runs over all temperature-dependent processes i with rate 
constants k,, where the temperature dependence of the rate 
constants is expressed in terms of the Arrhenius equation 
ki = Ar exp( — Ej/RT) [98]. Ai is the so-called pre-exponential 
factor and can, to a first approximation, be treated as temperature- 
independent. Eq. 1 1 applies to any kinetic model as long as the 
temperature dependence of the individual reactions are formulat- 
ed in terms of the Arrhenius law. 

The condition for temperature compensation is obtained by 
setting Eq. 1 1 to zero. Because in oscillatory systems the C[ 's have 
generally positive and negative values, there is a large set of 
balancing is, combinations which can lead to temperature 
compensation. The various combinations can be considered to 
arise by evolutionary selective processes acting on the activation 
energies [99]. Because the temperature homeostasis of circadian 
rhythms involves a compensatory mechanism [100], which needs 
to be distinguished from temperature-independence where all 
Cj^'s are zero, temperature compensation implies that there is a 
certain set of non-zero control coefficients with associated 
activation energies which (under ideal conditions) will satisfy the 
balancing condition Cj = 0 within a certain temperature range. 

The argument has been made that the balancing condition 
Cj = 0 should be non-robust and should therefore not match the 
many examples where mutations have no influence on the 
circadian period [101]. However, it should be noted that Eq. 11 
is model-independent and provides a general description how the 
period of an oscillator will depend on temperature in terms of the 
individual reactions defined by the A:,'s. Robustness, on the other 
hand, is a property of the actual oscillator model, where the 
number of zero C^'s can be taken as a measure for robustness. For 
the frequency controlled oscillators described earlier, there are 
certain regions in parameter space such as the shaded region in 
Fig. 7a, for which the oscillator's period is independent towards 
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